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Prologue 



Chiral Perturbation Theory (CHPT) is nothing but the low-energy limit of the 
Standard Model (SM) or, to be more precise, the effective quantum field theory 
describing hadronic interactions according to the SM, below the breaking scale of 
chiral symmetry (£ < A x ~ 1 GeV). This theory, founded a long time ago by the 
pioneering works of Weinberg, HD and Gasser and Leutwyler, @) is nowadays a rather 
mature subject. On one side several two-loop calculations have been performed in 
the purely mesonic sector, reaching, in some cases, a very high degree of precision. 
On the other side the original formulation has been successfully extended in several 
directions, including, for instance, heavy quark fields, bound-state dynamics, non- 
zero temperature effects, etc. 

A complete overview of the subject would be a tremendous task, and we 
could certainly not provide it in the four hours that we were given for these lectures. 
The interested reader is referred to some excellent reviews 3' @' 0) for a broader 
survey of the subject. The purpose of these lectures is to provide a basic introduction 
to CHPT. We will therefore restrict our attention to the simplest case, namely the 
mesonic sector, discussing in detail the determination of the effective Lagrangians 
and the calculation of a few specific quantities. 

The lectures are organized as follows: motivations, basic principles and 
the tree-level structure of CHPT are presented in the first lecture. Various aspects 
of loop calculations are discussed in the following two lectures. In particular, the 
renormalization at the one-loop level is presented in the second lecture, whereas 
the third lecture is devoted to studying the properties of unitarity and analyticity, 
and to showing how dispersion relations can be combined with the chiral expansion. 
Finally, the issue of non-leptonic weak interactions is introduced in the last lecture. 

1 Generalities and lowest-order Lagrangians 

1.1 Effective quantum field theories 

Within the Standard Model the interactions between quarks and gluons, ruled by 
Quantum Chromo Dynamics (QCD), are highly non-perturbative at energies below 
the breaking scale of chiral symmetry. This makes very difficult any description 
of the low-energy hadronic world in terms of partonic degrees of freedom. On the 
other hand, the spectrum of the theory is rather simple at low energies, containing 
only the octet of light pseudoscalar mesons: 7T, K and rj. Experimentally we also 
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know that, at very low energies, these pseudoscalar mesons interact weakly, both 
among themselves and with nucleons. It is then reasonable to expect that QCD can 
be treated perturbatively even at low energies, provided a suitable transformation 
of degrees of freedom is performed. This is exactly the goal of Chiral Perturbation 
Theory, where the pseudoscalar mesons are assumed to be fundamental degrees of 
freedom. 

Having an intrinsic energy limitation and being the low-energy limit of 
a more fundamental theory, CHPT is a typical example of effective quantum field 
theory (EQFT), a widely used tool in modern physics. The basic principle of any 
EQFT is that, in a given energy range, only few degrees of freedom are relevant and 
need to be described by dynamical fields. The remaining degrees of freedom of the 
more general theory can be integrated out, leading to effects that are encoded in the 
coefficients of appropriate local operators. 

Assuming this general point of view, all known quantum field theories, 
including the Standard Model, can be considered as effective. An important dis- 
tinction, however, is provided by the degree of renormalizability. In general the 
requirement of renormalizability, understood in the classical sense, is not manda- 
tory within an EQFT. H) Indeed if the theory is meant to be valid only for energies 
below a given cut-off A, and we perform an expansion of the physical amplitudes in 
powers of E/A, we can impose the weaker condition that, for any n > 0, the number 
of counterterms needed to regularize the amplitudes and contributing at 0[(E / A) n ] 
is finite. This condition is sufficient to ensure a predictive power to the theory for 
E < A. In other words, within an EQFT we only need renormalizability order by 
order in the energy expansion. The classical requirement of renormalizability, is a 
stronger constraint, which prevents a naive derivation of the intrinsic cut-off of the 
theory. The SM belongs to this restricted subclass of EQFTs without an intrinsic 
cut-off, whereas CHPT belongs to the more general case. 

1.2 Chiral symmetry 

Neglecting light-quark masses, the QCD Lagrangian can be written as 

£qcd = E qYU.-gs^GAq-^G^G^ + Oihewy quarks). (1) 

q=u,d,s \ / 

Apart from the SU(3)c local invariance, Cqcd nas a global invariance under the 
group SU(Nf)i x SU(Nf)fi x U(l)v X U(1)a, where Nf is the number of massless 
quarks (Nf = 3 in the case considered above). The U(l)y symmetry, which survives 
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also in the case of non-vanishing quark masses, is exactly conserved and its generator 
is the baryonic number. The U(1)a symmetry is explicitly broken at the quantum 
level by the Abelian anomaly. ® Finally G = SU {3)l x SU {3)r is the group of chiral 
transformations: 

( a \ 

^l,r — ► 9 L , R ^ L , R , where ip = I d I and g L R £ G . (2) 

If the operator ipif) has a non-vanishing expectation value in the ground 
state ((Ol^^lO) 7^ 0), or in the presence of a non-vanishing quark condensate, chiral 
symmetry is spontaneously broken. The subgroup that remains unbroken after the 
breaking of G is H = SU(3)v = SU(3)l + r, the famous SU(3) of the eightfold 
way, S) or, in the limit where only two quarks are kept massless, the SU(2) group 
of isospin transformations. 

The fundamental idea of CHPT is that, in the chiral limit (m u = irtd = 
m s = or, eventually, m u = m d = 0), the light pseudoscalar mesons are the 
Goldstone bosons generated by the spontaneous breaking of G into H (in the SU(3) 
case, m u = = m s = 0, the full octet of pseudoscalar mesons is identified with 
Goldstone fields, whereas in the SU(2) case, m u = = 0, only the three pions 
are). These light particles have the correct quantum numbers to be associated with 
the generators of G/H, as required by the Goldstone theorem. Moreover, since 
Goldstone fields can always be redefined so that they interact only through derivative 
couplings, 0) this hypothesis justifies the soft behaviour of pseudoscalar interactions 
at low energies. If pseudoscalar mesons were effectively Goldstone bosons, they 
would had been massless. This does not happen in the real world, owing to the 
light-quark mass terms, which explicitly break G. Since vn u ^^ s < A x , it is natural to 
expect that these breaking terms can be treated as small perturbations. The fact 
that pseudoscalar-meson masses are much smaller than the typical hadronic scale 
indicates that also this hypothesis is reasonable. Clearly this approximation works 
much better in the SU(2) case (Ml/M 2 p ~ 0.03) than in the SU(3) one (Ml/M 2 p ~ 
0.4). Summarizing, the two basic assumptions of CHPT are the following: 

1. In the chiral limit the SU(3)l x SU(3)r symmetry of the QCD Lagrangian is 
spontaneously broken into SU (3) l+r and the pseudoscalar meson fields can be 
identified with the corresponding Goldstone bosons. 

2. The mass terms of light quarks can be treated as small perturbations around 
the chiral limit. 
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According to these hypotheses, in order to describe the QCD interactions of pseu- 
doscalar mesons it is necessary to construct, in terms of Goldstone-boson fields, the 
most general Lagrangian invariant under G and add to it the explicit breaking terms 
(light-quark masses) that transform linearly under G. 0> The Lagrangian built in 
this way necessarily contains an infinite number of operators. Nevertheless, as antic- 
ipated, only a finite number of operators contribute at 0[(E/A x ) n ]). Therefore the 
theory can reach an arbitrary degree of precision for processes occurring at E < A x , 
provided a sufficient (but finite) number of couplings is fixed by experimental data. 

1.3 Non-linear realization of chiral symmetry 

Denoting by V{ the generators of H and by Ai the remaining generators of G, any 
element of G can be unambiguously decomposed as g = e^ Ai e ViVi . The Goldstone- 
boson fields are associated to the coordinates of the coset space G/H. To un- 
derstand how these transform under G we consider the action of a generic element 
g E G on u(£i) = e^ Ai : 

geG ge^ Ai = <V^ n ' . (3) 

The transformation — > provides a non-linear realization of the group 

G. 0) This realization is not linear since V^'s and A^s do not commute ([V^, Aj] ~ 
Ak); however, it becomes linear if restricted to the subgroup H: 

h = E H , h et* A * = U 0v V' A 'e-^l e^ Vi . (4) 



As shown by Callan, Coleman, Wess and Zumino, this non-linear realization 
provides the most general tool to construct operators that transform linearly under 
G (or invariant operators) in terms of the Goldstone-boson fields generated by the 
spontaneous breakdown of G into H. 

In the specific case of chiral symmetry an important additional information 
is provided by the automorphism induced on G by parity [P : (A iy Vi) — > (—Ai, VI)]. 
This implies that if gn : u(£i) — > u(^), then gi : u(— &) = — > u(ffi. We can 
therefore write 

«(&) 9 R u{ii}hr l (g^i) = h{g,£i)u{£^g~ l , 

u{0 gM&h- 1 (g,tt = h(g,ttu(Z i y g - 1 , (5) 

where h(g,£i) = e n (9&)v _ At ^his p 0m t we have all the ingredients to build generic 
operators transforming linearly under G, starting from their projection into H . For 
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instance, given a generic field \I/ transforming linearly under H as ^ e'^^e 'K^*, 
in the non-linear realization of G we find 

^^h(g,^h-\g,^) , (6) 

thus any product of the type (u, u*) x \I> x (m, u') transforms linearly under G: 

uW ^(wW)^ 1 , 

u t^ u _^ g L (u^u)g-\ 

m^m — ► g R (u^u)g~ , 

« f W g L (u^)g-\ (7) 

The above procedure can be generalized to the case of fields iff' belonging to different 
representations of H. Furthermore, starting from the derivative of u and we can 
define the following operators: 

Ufj, = i^d^u — ud^iv) = iwd^Uu' = , (8) 

rv = ±(u%u + ud lt rf) = -r] l . (9) 

It is easy to verify that both and the covariant derivative of 

y^ = d^-[Y^} , (10) 

transform according to Eq. (f|). Thus, similarly to Eq. (0), with appropriate con- 
tractions with u and we can construct operators transforming linearly under G 
also starting from and XJ ^ ■ Concerning the construction of terms invariant 

G — i 

under G, we note that given a generic operator Q transforming as Q — > giQdL 

G — i 

or Q — > gnQgn , an invariant term is obtained by the trace of Q in flavour space, 
denoted in the following by (Q). 

The choice of coordinates in the coset space G/H is not unique; however, 
in any given set of coordinates we can introduce a field u(£i) transforming as in 
Eq. (D.H) 

The freedom in the choice of coordinates implies that the parametriza- 
tion of u in terms of the pseudoscalar meson fields is not unique. In the following 
we shall adopt the exponential parametrization in the 3x3 flavour space, defined 
by 

u 2 = U = e^ F , 
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(ii) 



where rjs denotes the octet component of the r? meson. The parameter F appearing 
in Eq. ([TT]) is a dimensional constant (dim[F] = dim[$]) that, as we shall see in the 
following, can be related to the decay constant of pseudoscalar mesons. 

1.4 Lowest-order (strong) Lagrangian 

In the absence of external fields it is impossible to build non-trivial invariant oper- 
ators in terms of u and only, without their derivatives: it is necessary to have at 
least two derivatives acting on u or v) in order to build a non-trivial structure invari- 
ant under both chiral and Lorenz symmetries. If only two derivatives are considered 
this structure is unique: 

(«X> = (dpUPrf) . (12) 

Fixing the coupling constant of the operator (u^u^) so as to reproduce the correct 
kinetic term of spinless fields, leads to 



C 



(2) 



,(0) 



(13) 



This Lagrangian is the chiral realization of £>qcd a ^ the l° wes t order in the derivative 
expansion. 

~(2) 

Expanding C s in powers of $ leads to an infinite series of operators, whose 
couplings are all determined in terms of F: 

F 2 



(2) 



-{dJJd^) 



(14) 



~(2) 

From C s one can therefore determine the amplitude for any process of the type 
ix\i\2 ■ ■ ■ K n — > ^'1^2 ■ ■ ■ n 'm m the chiral limit. For example, considering the second 
term in Eq. ([14]) it is easy to show that 

(P+-P+) 2 



A{7T + 7l° -> 7T+7T ) 



[m,=0] 



F 2 



0(p 4 ) . 



(15) 



To parametrize the breaking terms induced by quark masses, and also to 
generate in a systematic way the Green functions of quark currents, it is convenient 



7 



to insert appropriate external sources both in the QCD Lagrangian and in its chiral 
realization. Following Gasser and Leutwyler, @) we introduce the sources 
and p, which transform as 



G 
G 
G 
G 



S + ip 
s — ip 

and we consider the Lagrangian 

Cqcd'- « - ^ - ^ (0) 

In this way we reach two interesting results: 
• The generating functional 



g R (s + ip)g~ 1 , 
g L (s-ip)g-\ 



(v, a, s, p) = C ( Q ' CD + ^7 M (^ + VflOV' - ^{s ~ ipTs)^ 



iZ(v,a,s,p) 



J VqVqVG e iJd 4 xC QCD (v,a,s,p) 



(16) 



;i7) 



is explicitly invariant under chiral transformations, but the explicit breaking 
of G can directly be obtained by calculating the Green functions, i.e. the 
functional derivatives of Z(v, a, s,p), at 



p = 



s = M g = diag(m M , m d ,m s ). 



(19) 



The global symmetry G can be promoted to a local one by modifying the 



transformation laws of and r M in 



r fJL = v fJL + 



Vfj. - a M 



— ► ^l^^ 1 + WLd^gl 1 . 



(20) 



Then the gauge fields of electroweak interactions can be automatically included 
as external fields by the substitution 

r 



a,. 



-eQA^ 
9 



g 



where 



2 cos 9w 
1 



Q cos(2^vi/) — 
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2 cos 9\y 



Q - 



6 
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2V2 



T + W+ + h.c. 








( 

000 
V 



(21) 



(22) 



and Vij denote the elements of the Cabibbo-Kobayashi-Maskawa matrix. 

As a consequence, Green functions for processes with external photons, Z or 
W bosons, can simply be obtained as functional derivatives of Z(v, a, s,p). 

The chiral realization of Z(v, a, s,p) at the lowest order in the derivative expansion 
is given by the classical action obtained from C s , after having included the external 
sources in a chiral invariant way. Concerning spin-1 sources, this can be achieved 
by means of the minimal substitution 

d„U -> DfJJ = d^U - ir^U + iUlp (23) 

or, in principle, introducing new operators written in terms of the tensors 

Fjf = dn v - d v F - i [F, l u ] , 

Fg = 9"r" -Vr 1 * -i[r^,r v ] . (24) 

From Eq. (^) it appears convenient to assign the same power counting to derivatives 
of u and to the sources a M and v M , so that D^U becomes a homogeneous term of 
first order in the derivative (or chiral) expansion: 



u ~ o( P { 



According to this assignment, the Lorentz-invariant terms containing the external 
tensors F^ V R are at least of 0(p 4 ) and do not appear at the lowest order. 

Regarding spin-0 sources, the most natural power-counting assignment is 
given by: 0) 

s, p ~ 0(p 2 ) . (26) 
As we shall see in the following, this choice is welcome since it implies M 2 ~ 0(p 2 ) 
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and leads to the Gell-Mann-Okubo mass formula. 

We are now able to write down the most general Lagrangian invariant under 
G, of order p 2 , which includes pseudoscalar mesons and external sources. This is 

4 2) = ^{Djjrrrt + x u ] + u x ] ) , (27) 

where 

X = 2B(s + ip). (28) 
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(2) 

C s is completely determined by chiral symmetry but for the couplings F and B, 
which have to be constrained from experimental data. These two couplings are in 
turn related to two fundamental order parameters of the spontaneous breakdown 
of chiral symmetry: the pion decay constant F n , defined by {0\ip'y' J "'y5ip\ 7l+ (p)) — 
i\[2F 1 ,p il , and the quark condensate (Q\tptp\0). Indeed, differentiating with respect 
to the external sources, we find 

5S {2) 

(0|^ 75 ^|vr + (p)) = (0\ \ 1T + (p)}=tV2Fp\ (29) 

A<?(2) 

(0|^|0) = -<0|-^-|0> = -i^S , (30) 

where is the classical action: 

S^ = fd A xCf . (31) 

It is important to stress that relations (p9|-|3~0D are exactly valid only in the chiral 
limit: in the real case (m q ^ 0) they are modified by 0(m 2 ) ~ 0(p 4 ) terms. 

The pion decay constant is experimentally known from the process 7r + — > 
fi + u: F n = 92.4 MeV. In principle one could determine F also from the kaon decay 
constant, defined analogously to F n and measured to be F K = 114 MeV. The differ- 
ence between F v and Fk is a typical 0(p 4 ) effect, which goes beyond lowest order. 
However, since 0(p A ) effects are expected to be larger in the case of Fk, the most 
natural determination of F at 0{p 2 ) is provided by F w . 

Contrary to the decay constants, the quark condensate is not directly re- 
lated to any physical observable. It is the product B x m q that can be constrained 
by means of experimental data. This appears in the quadratic terms of and is 
therefore related to pseudoscalar meson masses: 

Ml + = (m u + m d )B , 
M K + = (m u + m s )B , 

M K0 = {m d + m s )B . (32) 

The analogous equation for M 2 g contains no new parameter and gives rise to a 
consistency relation: 

3MJ = AM 2 K - Ml , (33) 

the famous Gell-Mann-Okubo mass formula, B ED well satisfied by experimental 
data under the assumption M v = M m . The validity of the Gell-Mann-Okubo rela- 
tion provides a significant a posteriori check that the 0(m 2 ) corrections to Eqs. 
are small and that the assignment s ~ 0(p 2 ) is consistent. 
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1.5 Quark- mass ratios and isospin breaking 

Similarly to the quark condensate, also the absolute value of light-quark masses 
cannot be determined within CHPT. Relations fl32|) provide, however, stringent 
constraints on quark mass ratios. 

In addition to 0{m 2 ) corrections, relations ( |3~2"D can be affected by elec- 
tromagnetic effects. At leading order in the chiral expansion, the latter can only 
depend on meson charges, and we can therefore write 

M 2 = (m u + m d )B + 0{m\) , 

M 2 K0 = (m d + m s )B + 0(m 2 q ) , 

M 2 + = (m u + m d )B + aA e _ m , + 0(m q ,am q ) , 

M 2 K+ = (m u + m s )B + «A c . m . + 0(m 2 q , am q ) . (34) 

Neglecting the small 0(m 2 , am q ) corrections and using the experimental values of 
pseudoscalar meson masses, one can extract from Eqs. ([34]) the following two ratios 

- {M «° ~ M " +) ~ {M " - M '* ] = 0.22 + 0.07 = 0.29 , (35) 



m d + m u M° 

m s (M K + - M,+) + (M K o - M w o 



25 . (36) 



m d + m u M° 

Interestingly, the three light-quark masses turn out to be rather different: Eq. (|35|) 
indicates a sizable violation (~ 30%) of isospin symmetry and Eq. (^) shows that 
SU(3) is not at all a good symmetry for quark masses. However, it is known that 
both symmetries, especially the isospin one, are usually well respected by strong 
interactions. For instance, considering the spectrum of vector mesons, it is found 
that isospin works at the (1-2)% level [(M w - M p )/M p ~ 1.5%] and 377(3) at the 
(10-20)% level [(M K * - M p )/M p ~ 16%]. 

The reason of this behaviour can be traced back to the smallness of quark 
masses with respect to the scale of chiral symmetry breaking. Indeed, although 
we cannot have a precise information about the absolute value of quark masses, by 
looking at the breaking of isospin and SU(3) symmetries in the spectrum of hadrons 
we can infer the following hierarchy: 

A x ~ 10 3 MeV , 
m s ~ 10 2 MeV , 
m d , m u ~ 10 1 MeV . (37) 
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The masses of all light hadrons, except for the would-be Goldstone bosons, remain 
different from zero in the chiral limit and are mainly determined by A x . Then isospin 
and SU(3) symmetries are accidental consequences of the fact that m Uj d, s *C A x . The 
observation that both these symmetries are well respected in the hadronic world is a 
good a posteriori check of the basic assumptions of CHPT discussed in Section |L~2} 

1.6 Lagrangian at order p 4 

The Lagrangian (E7p is only the first term of an infinite series. At the moment, in the 
purely strong sector the Lagrangian is known up to and including order p 6 . 0) Here 
we introduce the Lagrangian of order p 4 that was originally written down by Gasser 
and Leutwyler. @) In order to derive this Lagrangian, one needs to list all possible 
chiral-invariant terms of order p 4 . It is not strictly necessary, but very useful, to 
also reduce this list to a minimum, using all possible trace identities (which depend 
on the number of light flavours) and also the classical equations of motion derived 
from the Lagrangian (|27|), i n all possible ways. The interested reader is referred to 
the original articles, for more details on how this is done. Here we simply 



write down the Lagrangian in the case of three light flavours: 

£ 4 = L^D^D^U) 2 + L 2 (D^D V U){D^D»U) 

+L 3 (D^D»UD v U j D u U) + L 4 (D^D»U)( X j U + X U ] ) 
+L 5 (D^D»U( X j U + rf X )) + L 6 ( X j U + X U j ) 2 + L 7 ( X j U - X U ] ) 2 
+L 8 (x j U X j U + X U ] xU^) - iL^FgDpUDjfl + FjfDjflDJJ) 
+L 10 (U^F^UF Lfll/ ) + LniFR^Fg + F LfMV FD + L 12 (^x) , (38) 

where, for completeness, we have also written down the so-called "contact terms", 
those multiplying the constants L n and L 12 , which contain only external fields. We 



should add that at this order also the Lagrangian of Wess, Zumino and Witten 0) 
enters. This describes the effects due to the axial anomaly, which are therefore 
parity-violating. In the following we will not consider such effects. 



2 One-loop graphs: renormalization 

In the following two sections we will introduce the subject of loop calculations in 
CHPT. Technically, these loop calculations need no special introduction: any grad- 
uate student who has already made loop calculations (in QED for example) should 
be able to perform them also in CHPT also. On the other hand, at a conceptual 
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level, he/she may have doubts about the meaning of these loop calculations: how 
should one interpret the loops of pions, which contain contributions of pions of any 
virtuality, if one knows that at high energy scales these degrees of freedom are no 
longer relevant? The main aim of these two lectures is to answer these kinds of 
questions, and to illustrate the physical meaning of the loop contributions. 

2.1 The scalar form factor of the pion 

To illustrate the basic concepts of loop calculations in CHPT we will focus on one 
specific example: the scalar form factor of the pion, defined as 

(7r l ( Pl )^(p 2 )\m(uu + dd)\0) =: 5^T(t) , t = ( Pl + p 2 ) 2 , (39) 

where rh = (m s + m^jl. This matrix element is relevant to the decay h — > mi, 
which would have been the main decay mode for a light Higgs (of course this scenario 
is now experimentally excluded). The tree- level calculation of this matrix element 
is simple: in the Lagrangian (|27D the coupling between the scalar source and two 
pions, which can be read from Eq. (|28|) , does not involve derivatives and leads to 

T(t) = 2mB = M 2 + 0(p 4 ) . (40) 

This result, which we worked out from the Lagrangian, is actually a consequence of a 



general theorem, due to Feynman and Hellman. HID This states that the expectation 
value of the perturbation in an eigenstate of the total Hamiltonian determines the 
derivative of the energy level with respect to the strength of the perturbation: 

dM 2 

m— £ = (vrlrngglTr) = IYO) . (41) 
dm 

The value of the form factor at zero momentum transfer is fixed by this theorem, 
and a simple power counting implies that at leading order the scalar form factor 
is a constant - at order p 2 the theorem completely fixes the form factor. On the 
other hand no general principle forbids a dependence of the form factor on t (to the 
contrary, they imply it, as we will see), and to generate this we necessarily have to 
go beyond leading order. 

2.2 SU(2) Lagrangian at order p 4 

Before starting the loop calculation, let us have a look at what happens to the 
form factor once we include tree-level contributions from higher-order terms in the 



Lagrangian. The Lagrangian at order p has been discussed in Sect. 1.6 in the case 
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of three light flavours. For the case at hand, the role of kaons and etas is marginal, 
if we stick to the very low-energy region. It is more convenient to use a simpler 
Lagrangian, by expanding around m u = m d = and keeping m s at its physical 
value. In this case the chiral symmetry is SU(2)l x SU(2)r. The Lagrangian of 
order p 2 remains unchanged - the only change is that the field U, and its logarithm 
<3> are now 2x2 matrices. At order p 4 the Lagrangian is simpler, because we can 
use more trace identities to reduce the number of independent terms. For two light 
flavours this is a sum of seven terms^ 

Cf = Y.kP l + ... , (42) 



where 



=i 



Pe = \(K,u u }fn , P 7 = -^( X -} 2 , (43) 
and we have used the compact notation: 

X± = u] XU ] ± u\ ] u 

= uFl v u ] ±u ] F^u . (44) 

Only two of these seven terms contribute to the scalar form factor: the terms pro- 
portional to £3 and Z4. Their contribution reads: 

M 2 r 1 

T [hM (t) = — [-4M% + th] ; (45) 

the calculation is recommended as an easy exercise, as is the calculation of the 
contribution of Z3 to the pion mass. Once these two calculations are completed, one 
can then check that the Feynman-Hellman theorem is respected also in this case. 

The result in Eq. (f45[), a tree- level calculation with the next-to-leading 
order Lagrangian, is merely a statement about how the symmetry constrains this 
particular matrix element: at this order the scalar form factor can have at most a 
term linear in t. No symmetry relation exists between the constant term and the 

2 The ellipsis is for the so-called contact terms, i.e. those that depend only on the external 
sources. 
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Figure 1: One- loop diagram for the scalar form factor of the pion. The double line 
stands for the scalar source, whereas single lines for pions. 

coefficient of the linear term, hence we have two different constants. The constant 
term is related to the derivative of the pion mass with respect to the strength of 
the symmetry-breaking term in the Lagrangian, whereas the coefficient of the linear 
term is related to the correction to the pion decay constant (again, it is a very good 
exercise to calculate the latter with the Lagrangian in Eq. (|42|)). 

2.3 Loop graphs 

If we neglect the tadpole graphs and those renormalizing the external legs, there is 
basically only one graph for this process, the one shown in Fig. [I]. Its structure is 

r dH {M 2 ,p 2 ,p-l,l 2 } r dH 1 

J (2tt) 4 (/ 2 -M 2 )(( p -/) 2 -M 2 ) ^ X J (2tt) 4 (/ 2 -M 2 ) 

r dH 1 
+ V J (2tt) 4 (/ 2 -M 2 )((p-/) 2 -M 2 ) 
= xT(M 2 ) + yJ(t) , (46) 

where p = p\ + P2- We have indicated, in the first integral, all the terms that 
can appear in the numerator, and, after the arrow, the two possible structures to 
which the various terms can be reduced. The momenta and masses in the numerator 
come from the four-pion vertex on the right-hand side of the diagram. The power 
counting for this integral shows that it represents a correction of order p 2 to the 
leading-order term (the integration measure, which is of order p 4 , is compensated 
by the two propagators). This is true in general: one- loop graphs constructed from 
the Lagrangian (^7|) are always a correction of next-to-leading order in the chiral 
expansion. For example, it is easy to see that no matter how many vertices we add 
on the internal lines in the graph in Fig. |l], the graph remains of order p 2 : if we add 
a vertex of order p 2 we also get an extra propagator, which compensates it. The 
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general theorem was first proved by Weinberg. E) We will come back to this point 
later. 

If we expand the tadpole integral T(M 2 ) and the loop integral J(t) in a 
Taylor series in their respective arguments: 

T(M 2 ) = a + bM 2 + f(M 2 ) , 

J{t) = J(0) + J(t) , (47) 

only the first terms in the expansion are divergent, whereas both T(M 2 ) and J(t) 
are finite - this is easily seen by taking a sufficient number of derivatives on the loop 
integrals in Eq. (fR]) . It is left as an exercise to the reader to show that the one-loop 
diagrams that we have neglected can only produce terms like T(M 2 ). Therefore 
the contribution of the loop diagrams to the scalar form factor has the following 
structure: 

M 2 r - i 

r(t) ~ —» x x bM 2 + x 2 t J(0) + xif(M 2 ) + x 2 J(t) . (48) 

The divergent part of the loop diagrams has exactly the same structure as the 
counterterm contribution calculated above: to remove it we simply need to define 
the counterterms properly (in this case the constants Z 3 and I4). 

The principles stated in the first lecture that guided the construction of 
the effective Lagrangian only appealed to symmetry arguments, and therefore al- 
lowed for an infinite number of terms. Once this is accepted, there is no problem 
of principle in carrying through the renormalization program: as anticipated, the 
difference between renormalizable and non-renormalizable theories is, in a sense, a 
technical detail. 

2.4 Chiral logarithms 

If we expand the form factor in a Taylor series in t, we can write it in the following 
form: 



r(f) = r(o) 



l + \(r 2 y s t + 0(t 2 ) 



(49) 



The coefficient of the linear term, properly normalized, is called the scalar radius of 
the pion. Its size is a way to represent the spatial extension of the pion when probed 
with a scalar source. We have stated above that the coupling constant that appears 
in this quantity, l±, also determines the first correction of the pion decay constant 
around the chiral limit. There is another piece of information on the scalar radius, 
which we can already gather from the simple sketch of the loop calculation given 
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above: 

^-^-Iwtw^mt^' (50) 

namely that the scalar radius contains an infrared divergence. In the chiral limit this 
quantity diverges. This divergence should not be removed and does not represent a 
problem because it has a physical meaning, in the sense that when the pion becomes 
massless the cloud of pions surrounding any hadron (and therefore also the pion 
itself) extends to an infinite range. A quantity that measures the spatial extension 
of this cloud should indeed become infinite in the chiral limit. Notice that the scalar 
form factor is finite and remains finite also in the chiral limit - it is only the first 
derivative in t, calculated at t — that diverges when the pion mass goes to zero. 

These infrared divergences are present everywhere in pion physics, and 
in many cases they are among the most important physical effects (less so in the 
case of kaons). Their relevance was first pointed out by Li and Pagels. IHP The 
effective Lagrangian method provides a systematic way to calculate these effects. 
As we have seen in the above example, they arise from the infrared region in the 
loop integrals - precisely the region where we should fully trust the vertices of our 
effective Lagrangian. 

One may be less at ease with the ultraviolet region of the loop integrals: 
there indeed has no justification the use of the effective Lagrangian. On the other 
hand, through the process of renormalization, that part of the integrals is completely 
removed and substituted with unknown constants, the counterterms. As is some- 
times said, these parametrize our ignorance of the physics that lies above the range 
of applicability of the effective Lagrangian. In fact, we are not completely ignorant 
about the physics of strong interactions in the GeV range and above, and as we will 
see in the following lecture, it is not difficult to get a good estimate of the numerical 
value of the counterterms. 

Until now the only part of the loop integrals that we have not analysed is 
the finite, analytically non-trivial part of the loop integral, the function J(t). This 
will be treated in full detail in the following lecture. Before doing that, however, we 
want to consider the problem of renormalization from a more general point of view. 

2.5 Renormalization and chiral symmetry 

The reader who has done the exercise of calculating the divergent part of the scalar 
form factor to one loop will have seen that, to renormalize it, the two counterterms 
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Z3 and h have to be denned in the following way 

h = l r M-\\ , h = l r M + 2\ , (51) 

where A is divergent and dependent on the regularization method (in dimensional 
regularization, for example, it is defined by A = (cfi) d ~ A / (d — 4), with c an arbitrary 
constant that defines the regularization scheme). 

If we look at the definition of the operators in front of 1% and 1^ we see that 
(as usual) they contain an infinite number of pion fields. For example they both 
contain a term with one scalar source and 6 pion fields. Does it mean that if we had 
calculated that matrix element to one loop we would have found the same divergent 
part for l 3 and Z 4 ? Or, to put it differently, that once we calculate the divergent 
part of the scalar form factor then we know the divergent part of all other matrix 
elements with a higher number of pion fields^? 

The answer is yes. Chiral symmetry puts a strong constraint on the diver- 
gences: they have to be chiral-invariant terms. This conclusion is originally due to 
Weinberg, HP on the basis of a highly plausible, but still heuristic argument. It is 
now put on a solid basis by the work of Leutwyler: he proved that to calculate 



hadronic Green functions with an effective Lagrangian, such that they respect the 
Ward identities implied by the chiral symmetry, one necessarily has to start from 
a chiral-invariant effective Lagrangian. The non-trivial part in this statement is 
that it takes into account also quantum effects: anomalous symmetries show that 
it is not always true that a symmetry that exists at the classical level survives the 
quantum corrections - or vice versa, that to have a symmetrical quantum theory 
one necessarily has to start from a symmetrical classical Lagrangian. 

The theorem applies also to the divergent part of the quantum corrections: 
they have to be chiral-invariant. A general theorem of quantum field theory states 
that the divergent part of a loop graph is a polynomial in the external masses and 
momenta. These two theorems lead to the conclusion that the divergences, order 
by order, can be reabsorbed by the chiral-invariant counterterms. 

While these general theorems are very important, performing a direct cal- 
culation of the divergences and the subsequent renormalization in an explicitly 
chiral-invariant form, is probably much more instructive. The rest of this lecture 
will be dedicated to an illustration of the tools necessary for this calculation and to 
a sketch of the procedure. 



3 Strictly speaking, these questions make full sense only if we suppose that these two are the 
only counterterms at order p 4 . 
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2.6 Background field method and heat kernel 

To introduce these techniques, it is useful to consider a simpler theory than CHPT: 
a purely scalar O(N) 4 theory. Its Lagrangian is the following: 



71=1 

A, = \ (<9^W - MW) - f (0>*) 2 - , (52) 

where £o is the classical Lagrangian and Ci, i > 1, the series of counterterms needed 
to renormalize the theory, and a summation over repeated indices is implied. So as 
to make the situation completely analogous to CHPT, where we have the external 
sources v, a, s and p, we have also introduced external sources /' coupled to 
the fields 0\ A path integral constructed with this Lagrangian is a function of the 
external sources By taking the appropriate functional derivatives with respect 
to the sources we can extract from the path integral all possible Green functions 
of the fields 0\ The logarithm of the path integral is usually called the generating 
functional Z{f}: 



e iZ{f}/h 



= Nj [d<P]e iS l h , S = JdxC , 

= z + nz 1 + h 2 z 2 + o(h 3 ) , z{o} = o . (53) 



To calculate the classical part of the generating functional Z and its quantum 
corrections Z$, it is useful to expand the field around the solution of the classical 
equations of motion, 0: 

= (M 2 + D)0* + gpfi + P = • (54) 

We can then shift the integration variable: 

^ = ^ + f , [d4>] = [dt] , i-Oih 1 ' 2 ) , (55) 

and get, for the path integral: 

e iZ{f}/h = Ne ts/n J ^ exp || y dx i< '/;,/ ' + 0(£ 



(56) 



where 



Dij = —OSij + (Tij , 



{M 2 + gj> 2 )8 ij -2gP4? . (57) 
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The separation of the field into its classical part and its quantum fluctuations £ 
makes the calculation of the Taylor series of Z in H more transparent. Indeed one 
immediately sees that the classical action S n of the Lagrangian C n , evaluated at 
the solution of the classical equation of motion contributes only to the term Z n of 
the generating functional. The path integral over the quantum fluctuations yields 
higher-order quantum corrections. For example Z\ receives a contribution from the 
integral over the quadratic term in £ , which was explicitly given in Eq. (p6|). The 
path integral of the exponential of a quadratic term is known, and can be given in 
closed form - it is the determinant of the differential operator Dy: 



Z\ = J dx 



i 



-ln^ZVl+A 



(5? 



where D = D\ f=0 = -(a+M 2 )5 tj . To acquire familiarity with the formal expression 
(|58|) it is useful to expand the operator D around / = 0: Dy = D 0i j + rjij, and 
correspondingly expand the logarithm around 1: 

In (DyDoy 1 ) = %-A>y - -VijDojkVkiDou 1 + ■■■ ■ (59) 

It should now be clear that this logarithm simply represents the sum of all one-loop 
diagrams with any number of four-0 vertex insertions, as shown in Fig. [2|. The 
reader with some experience of loop calculation will have immediately realized that 
only the first two terms in the expansion of the logarithm generate divergences: all 
the others are finite loop integrals. Indeed the divergent part of Z\ can be calculated 
by working out two simple loop diagrams. 

The situation would not be as simple in the case of CHPT: there, the pres- 
ence of derivative couplings complicates the situation significantly. The calculation 
can nonetheless be done in close form with the help of the heat kernel method. 
Without giving many details let us see how it works in the case of the O(N) theory. 

A differential operator of the form: 

D 2 = -d 2 + a , ^ = <9 m + 7m , (60) 

where a and 7 are matrices in general, is called an elliptic differential operator, 
and has the same form as the differential operator that appears in the equations 
describing the diffusion of the heat (which explains the origin of the name). The 
divergent part of the logarithm of such operators is known in terms of the operators 
a and 7^, without knowing their explicit form. It is given by: 



Jdxln(D tJ D ^) = ^^^ Jdx 



o 



+ ... . (61) 
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Figure 2: Graphical representation of the one- loop contributions to the generat- 
ing functional. The logarithm of the differential operator, represented as the loop 
without external legs, is the sum of all diagrams with any number of external legs. 

Details on how this result is obtained can be found in Ref. E3) . In the case of the 
O(N) theory that we are considering, 7^ = and we simply have to calculate the 
square of a: 

a 2 = 2(N + 2)gM 2 (j) 2 + (N + 8)g 2 (f) 4 . (62) 

The calculation of the divergences is complete. The reader is urged to verify this 
result by an explicit calculation of the two divergent diagrams in Fig. |2|. 

2.7 Calculation of the one-loop divergences in CHPT 

With the tools introduced in the preceding section we have reduced the problem 
of the calculation of one-loop divergences to a simple algebraic exercise. Even in 
the case of CHPT. The only step that requires some care is the choice of the field 
£, representing the fluctuations around the classical solution. Indeed this can be 
done in many different ways, but only a few of them will simplify and make the 
intermediate steps of the calculation transparent. One of the most convenient choices 
is the following: 0) 

U = e^ F , = + £ , 

U = ue ii/F u . (63) 

As we have seen in the first lecture, the transformation properties of the exponential 
of the £ field are particularly well suited to check the chiral invariance properties 
during the various steps of the calculation. 

The expansion of S2 = / d^xCg around the classical solution reads 

f (2) f F 2 

J dxC s = J dx—(u fl u fl + x+) 
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J dxcf + J dx l -c^e + o(e) , (64) 



where 



A = -d 2 + a , 

d/i, ij d^dij -\- 7^, ij , 
1 
2 



= --([A i; Ajjr^) , 



r, 



Oij = —^(Wti,M[^j,u^) + {Xi,Xj}x+) • (65) 

It is now a simple algebraic exercise to calculate the divergent part of the 
generating functional to one loop (e = d — 4): 

-^logdetA = ^ljd d x[^(([u,,u u ]-2if + , u ) 2 ) + ^((uX + X + ) 2 ) 

This expression is valid for a generic number of light flavours Nf, and is explicitly 
chiral-invariant, as it should be. To obtain the result in the interesting physical cases 
Nf = 2, 3 is not only necessary to trivially substitute Nf with the numerical value 
of interest, but also to reduce all the chiral-invariant terms appearing in Eq. ( |5^D 
to a minimal set. This step involves the use of trace identities and of the equations 
of motion. Without giving the details of this final step, let us conclude this lecture 
by presenting the results in the case of Nf = 2. If we define the counterterms as 

l i = r i ( J JL)+'y i X(jJL) , (67) 



their coefficients for the divergent part are: 

12 1 

7i g j 72 g j 73 g 



11 

74 = 2, 75 = --, 7e = --, 77 = . (68) 

6 6 
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3 One-loop graphs: analyticity and unitarity 

According to the Watson theorem, 0) above threshold but below the inelasticity 
threshold, the phase of the scalar form factor is equal to the S"-wave tttt phase shift 
with isospin 1 = 0, 5®. As is well known, this theorem is a consequence of unitarity: 

Imf(t) = a(t)T(t)t° *(t) = T(t)e- l5 °sm5° 

(69) 

= |r(*)|sin5 , (70) 

where a(t) = [1 - AM 2 /t} 1/2 and F(t) = T(t)/r(0), and tg is the I = S-wave of 
7T7T scattering. The unitarity relation (|70|) shows that the leading-order expression 
T(t) = 1 cannot be the whole story: if we want an accurate description of the form 
factor away from t = 0, we need to include higher orders and, in particular, loops - 
imaginary parts can only be generated by loop graphs. 

Notice that since at leading order T(t) is 0(1), and the phase Sq is 0(p 2 ), 
the 0(p 2 ) imaginary part (which is a next-to-leading order correction) is completely 
fixed by leading-order quantities: 

lmT^(t) = 6f\t) = a(t) 2 -^ . (71) 

The use of the effective Lagrangian method to calculate the form factor guarantees 
that this relation is satisfied. The complete expression for the one-loop scalar form 
factor reads as follows: 

t If — M 2 - 



where J(t) is the subtracted one-loop integral (f47P. Its explicit expression reads: 



At) 



16vr 2 



/ x , cr(t) - 1 

a(t) In -H + 2 

1 ; a(t) + 1 



(73) 



The reader can now easily verify that the imaginary part of the form factor at this 
order indeed satisfies (|7Tr). 



3.1 Dispersion relation for the scalar form factor 

An analytic function must be real on the real axis: the scalar form factor is non- 
analytic from threshold (4M 2 ) up to infinity. On the basis of very general arguments, 



23 



which mainly use the causality principle, one can prove that as a function of t, 
the scalar form factor must be analytic everywhere else (see Ref. E3) for a general 
discussion of the analyticity properties of amplitudes and Green functions). Since 
it is analytic everywhere else, the non-analyticity of the form factor on the real axis 
can be further characterized, and described as a discontinuity: 

T(t + is) = T*(t- ie) = |f (t)\e< {t) . (74) 

Given these analytic properties, we can write the following dispersion re- 
lation: 

t 2 f°° dt'lmT(t') . . 

r(t =1 + W + - /. — — V , (75) 



7T J4M| t' 2 t' — t 

where, for later convenience, we have subtracted the dispersive integral twice - we 
will come back to the issue of how many subtractions are necessary for the dispersive 
integral to converge. The dispersion relation shows that, if we know the subtraction 
constants (in this case only one, b) and the imaginary part on the real axis, we can 
reconstruct the scalar form factor everywhere on the complex plane. 

It is no surprise that any perturbative calculation in a quantum field theory 
produces amplitudes and Green functions with the correct analytic properties. Using 
an effective field theory makes no difference: the form factor calculated to one 
loop in CHPT must have the correct analytic properties, and must satisfy (at the 
perturbative level) the dispersion relation ( [T5|) . To convince ourselves that this is 
actually the case, let us first apply the chiral counting to the dispersion relation: 

f(o) = i 

b ~ 0(1) 

Imf(t') = 0{p 2 ) . (76) 

As we have seen the 0(p 2 ) imaginary part is fully fixed by leading-order quantities, 
([H|), and apart from an unconstrained polynomial term, the real part must be given 
by the dispersive integral over this known imaginary part. We leave it as an exercise 
to prove that this is true. For this it is useful to know that: 

j (t) « r *£<o . (77) 

v ; 16tt 2 Umi f t'-t y J 

3.2 High-energy contributions to the dispersive integrals 

In the previous section we showed that the renormalization procedure removes the 
contributions to the loop integrals where the momentum squared of the pions is 
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large. This was reassuring because we cannot hope that our effective Lagrangian 
describes highly virtual pions well. In the present section we are dealing with the 
contribution to the loop integrals from real pions: the dispersive integrals. As we 
have seen above, these extend all the way up to infinity, as required by analyticity. 
In the perturbative expansion that we are considering, the imaginary part of the 
form factor, which is in the integrand, is evaluated only to leading non-trivial order. 
This description of the imaginary part can be reasonably accurate only in the low 
energy region: still, in the integral, it is used all the way up to infinity. How can we 
trust the dispersive integral? 

In fact, we don't. At least not for the contribution that comes from the 
high-energy region. Suppose we decide to remove the part of the dispersive integral 
from A = 1 GeV to infinity. We should then subtract from the form factor a term 
like: 

The Taylor expansion inside the integral can be safely performed because the CHPT 
calculation of the form factor is valid only for i < A 2 . In the chiral language, this 
part of the dispersive integral can be represented as a polynomial series starting at 
order p 4 , i.e. at an order which is beyond the accuracy at which we are presently 
working. This shows that worrying about these contributions to the dispersive 
integrals is futile: the only sensible way to improve the evaluation of the dispersive 
integral is to go one order higher in the chiral expansion. This would automatically 
give a representation of the form factor that contains the dispersive integral over 
the imaginary part correct up to order p 4 . 

For those who are interested in the numerics, the first term in the Taylor 
expansion of the integral (|78|) is equal to: 

t 2 r^sf\f) 



7T J A* t' 3 



0.7 



t(GeV 2 ) 2 , (79) 



which means, for t = (0.5 GeV) 2 (which is about the upper limit of validity of the 
chiral expansion), a 4% correction to the leading-order result. Also numerically 
everything is well under control. 



3.3 Estimate of the low-energy constants 

The chiral representation to next-to-leading order satisfies a dispersion relation with 
two subtractions. The number of subtractions in this case is dictated by the be- 
haviour at infinity of the leading-order chiral phase 5®. The latter, however, has 
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nothing to do with the physical reality: it can be shown that, with form factors, one 
subtraction is already sufficient for the dispersive integral to converge. To discuss 
the physical implications of this, it is useful to consider the vector form factor of the 
pion, which is defined as follows: 



^ l (pi)\qr k l,qW(p 2 ))=ze^F v (t) , q= [ U . . (80) 




It satisfies the following dispersion relation: 

Fv(t) = i + 1 r ^Mm . (81) 

This dispersion relation implies that the charge radius (the analogue of the scalar 
radius, for the vector form factor) is given by a sum rule: 

r 2 Y v = - / dt' -P } . 82 

6 7T J AM* t 

If one knows the imaginary part of the vector form factor one can calculate the 
charge radius. 

Within the chiral representation the sum rule (|52p does not make sense as 
it stands, because the integral on the right-hand side does not converge. In CHPT, 
the charge radius is given in terms of one of the low-energy constants: 



Vr - -- 



S3) 



and, in this language, the sum rule ( |52"D can be read as a way to determine 
Notice that the fact that this constant appears at all in the chiral representation is 
a consequence of the non- convergence of the chiral dispersive integral: that integral 
would receive a sizeable contribution from the high-energy region, and there the 
chiral representation obviously fails, not even allowing the calculation to be made. 

Having the sum rule at hand, we can attempt to estimate this particular 
low-energy constant, from an estimate of the behaviour of the imaginary part of the 
form factor. What do we know about it? Actually, this has been measured in various 
ways, and a lot of experimental information is available on this quantity. But rather 
than collecting all the available information to make an accurate evaluation of the 
sum rule, we want to make a simple, but instructive exercise. Since in this channel 
(two pions in the P wave and isospin / = 1) the p resonance contributes, we can try 
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Table 1: Contributions of the resonances V, A, S, S\ and r/i to the constants L\ 
in units of 10~ 3 . 0) The phenomenological values in the second column are either 
from Ref. 0> @) (entries 4, 6, 9 and 10) or from Ref. (all the remaining ones). 



to give a rough estimate of its contribution to the sum rule. For this it is sufficient 
to use the narrow-width approximation: 

lmF v (t) ~ Ti5{t-M 2 p ) , (84) 

which, when inserted in the integral, gives 

lirTv = §^ , (85) 
where the two coupling constants 

F p =144MeV, / = 69MeV (86) 

determine the strength of the coupling of the p resonance to the external vector field 
and to the pions. 

These numbers give 

l r 6 (M p ) ~ -13.3 x 10" 3 , (87) 
whereas if we extract it from the measured value of the charge radius, we obtain: 

l r 6 (M p ) = -(13.5 ± 2.5) x 10~ 3 , (88) 
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in very good agreement with fl87|) . This is not a big surprise: we are simply com- 
paring the experimental determination of the charge radius with the contribution 
that the p gives to it, and it is well known that the latter is largely dominant. It 
is however illuminating on the physical meaning of the low-energy constants in the 
chiral expansion. Physically, these constants must encode the effects of the physics 
which occurs above the Goldstone-boson scale, and which does not appear explicitly 
in the effective Lagrangian. Intuitively one would expect that the resonances are 
the most important "high-energy" phenomenon that is neglected here. A simple 
estimate of their role should give account of at least the order of magnitude of these 
constants. 

Actually a careful analysis of the contributions of all the lowest-lying res- 
onances to the low-energy constants shows that they practically saturate the exper- 
imentally determined values. 0) This is clearly seen from Table [I]. 

3.4 Exact solution of the dispersion relation 

If we assume that only the tctt channel is open all the way up to infinity, then the 
phase of the form factor on the cut is everywhere the tttt phase shift <5q. In this 
approximation, we can pose the following mathematical problem: suppose that 

1. F(t) is an analytic function of t in the whole complex plane, with the exception 
of a cut for AMI < t < oo; 

2. approaching the real axis from above, e~ lS ^F(t) is real on the real axis, where 
5(t) is a known function. 

Can F(t) be determined in general? The solution to this problem is due to Omnes, H3) 
and reads as follows: 



( t r 

F(t) = P(t)Q(t) = P(t) exp{- \ 

[ 7T J 4 



t /•« dt' 5(f) 



4Ml f t'-t 
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where P(t) is a generic polynomial, and Q(t) is called the Omnes function. P(t) can 
only be constrained by the behaviour of the function F(t) at infinity: if we know the 
behaviour of the phase 5, and also that of F(t), we can fix at least the degree of the 
polynomial P(t). Then, depending on the degree of this polynomial, we will need a 
number of inputs to fix its coefficients. Whatever the degree of the polynomial, the 
relevant conclusion here is that if we know the phase, we only need to determine a 
few constants to fix the function F(t) completely. 
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Can we use this information and improve the chiral representation of the 
form factor combining it with the Omnes representation? 

In the chiral representation at order p 4 we have seen that there appears 
the Omnes function expanded to next-to-leading order: 

rg) pT (t) = 1 + b®t + A«(f) , A^\t) = f f ^ . (90) 

7T JAM% V t — t 

Indeed the above expression can be seen as a chiral expansion of the Omnes solution 



r(t) = (l + bt)Q(t) ~ (l + [6 (0) + 0(M 2 )]t) [l + A^ 2 \t) + 0(p 4 

= l + b {0 k + A (2) (t) + 0(p 4 ) , (91) 

and we could try to reconstruct the full Omnes solution from its expansion. An 
obvious improvement would be, for example: 

T Q (t) = (l + b^t)e A(2) ^ , (92) 

where we have simply exponentiated the dispersive integral over the 7iir phase; or we 
could use the chiral phase at next-to-leading order ID in the exponentiated dispersive 
integral: 

I^(t) = (l + &^)e A(4) «; (93) 

or, if we have a good phenomenological representation for the phase, we could use 
this to calculate the dispersive integral: 

fr(t) = (1 + fr(2) t ) e A phys « _ ( 94 ) 

Also, in calculating the dispersive integral, we could in principle choose whether to 
use a cut-off or to extend the integration up to infinity. However, given the definition 
of the Omnes function (|89|) , the number of subtractions is not sufficient to make the 
integral convergent, even for the p 2 phase: with a once-subtracted Omnes function 
using a cut-off is mandatory. 

Finally, we have to mention another possible degree of freedom in combin- 
ing the chiral and the Omnes representation: the choice of the subtraction point. In 
the present case there seems to be no choice, because we know that the form factor 
at t — must be equal to 1. It is then natural to have both the polynomial and 
the Omnes function equal to 1 at t — 0. This choice fixes the subtraction point of 
the Omnes function. If the information on the value of the form factor at zero were 
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Re(F*(s)/Fj;0)) 




Figure 3: Comparison of the physical scalar form factor (solid line) with various 
"Modified Omnes representations" (MOR). The dashed and dot-dashed lines show 
the difference between the use of the phase to one and two loops, whereas the long- 
dashed line shows that by tuning some of the low-energy constants, the MOR can 
actually reproduce the form factor in a wide range of energy. 



missing, we could as well have chosen a different subtraction point of order M%, 
without changing the chiral counting of the final result. 

There are various ways to combine the chiral and the Omnes representa- 
tions. In all cases we improve the chiral representation since we are able to sum up 
higher-order terms that are fixed by unitarity. On the other hand the difference be- 
tween one way and the other to implement this improvement is only one order higher 
than the chiral representation we started with. We can claim a real improvement 
only if we can show that this arbitrariness can be constrained by physical arguments, 
and that its effect is numerically small. For example, If we had a good phenomeno- 
logical representation for the phase, with very small uncertainties, it would certainly 
be preferable to use the representation f 7 (t) rather than T a (t), or T^(t). We have 
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already commented on the choice of the subtraction point: that is another example 
of an important information that helps in reducing the arbitrariness of this unita- 
rization procedure. In the present case the only degree of freedom is in the choice 
of the linear term of the polynomial: Eq. fl94"|) shows that combining the chiral and 
the Omnes representations amounts to fixing this subtraction constant with CHPT. 

In the recent literature there are various examples of the application of such 
procedures: one combines the solution of a dispersion relation (be it exact, as in the 
case of the form factor, or numerical, in other cases) with the chiral representation, 
and uses the latter only to determine the subtraction constant. The advantage of 
such a procedure is clear: the dispersive integrals, which extend all the way up to 
infinity, are certainly dominated by the low-energy region, but also receive sizeable 
contributions from the intermediate-energy region, between 0.5 and 1 GeV. Using 
the chiral representation to evaluate the contribution from that region is certainly 
not the best one can do. If one has better information on imaginary parts, or 
phases, in that region, one should definitely use it. As we have discussed here, this 
procedure is not free from arbitrariness, and all care should be used to constrain 
this to a minimum. The interested reader is referred to the original articles for 
a more detailed treatment of the case of the pion form factors. A review of the 
marriage between CHPT and dispersion relations can be found in Ref. 

3.5 Concluding remarks on lectures 2 and 3 

In these two lectures we have analysed in detail the various contributions entering 
the next-to-leading order calculation of a typical amplitude in CHPT. To make the 
analysis concrete we have considered the scalar form factor. 
We have seen that: 

1. The renormalization procedure is straightforward - indeed not very different 
from the case of a renormalizable field theory. 

2. The loop integrals generate infrared divergences that are physical: the effec- 
tive-Lagrangian method is a systematic way to calculate them. 

3. The finite, analytically non-trivial part of the loop integrals is dictated by the 
properties of unitarity and analyticity. As usual in quantum field theory these 
properties are built in, and automatically respected in loop calculations. 

4. The finite part of the counterterms can be expressed in terms of sum rules 
and in most cases these are saturated by the lowest-lying resonances. The 



31 



low-energy constants have a very clear physical meaning, as they embody the 
effect of the high-energy degrees of freedom that are not explicitly considered 
in the effective Lagrangian framework 



4 Non-leptonic weak interactions 

4.1 Partonic \AS\ = 1 effective Hamiltonian 

As discussed in the first lecture, the chiral realization of Z(v,a,s,p) allows us to 
calculate not only transition amplitudes of pure QCD, but also Green functions of 
weak and electromagnetic processes with external gauge fields, such as semileptonic 
kaon and pion decays. Z(v, a, s,p) is not sufficient, however, to describe non-leptonic 
weak transitions, such as K —>■ 2tt and K — > 3n decays. In this kind of processes 
the W boson is coupled to two quark currents and cannot be treated as an ex- 
ternal field. The strong-interaction effects that renormalize the non-leptonic weak 
transition cannot be trivially factorized and should be evaluated up to distances of 



on a twofold EQFT approach. The first step, performed within perturbative QCD, 
is the construction of the partonic | AS\ = 1 effective Hamiltonian, or the expansion 
of the non-local product of weak currents into a series of local partonic operators 



review). This first step, based on Wilson's Operator Product Expansion, allows 
us to encode in appropriate coefficients the sizeable QCD corrections that renor- 
malize the weak interaction at short distances, from Myy down to the lowest scale 
where perturbative QCD can still be applied, considerably simplifying the original 
problem. 

Denoting by J^{x) the charged weak current and by D(^(x, M w ) the W 
propagator in spatial coordinates, the partonic |Aif>| = 1 effective Hamiltonian at 



0(1/M W ). 



The simplest strategy to describe non-leptonic weak interactions is based 




.ensive 



0(G F ), 

n[ A f f =l = -^v* s v ud Ci(ji)Qi(ji) + h.c. , 



(95) 



can formally be defined by the following equation 



T(J -> F) 



1 



d 4 xD^(x, M W )(F\T (jj.Oc) Jj(0)) |/> + 0{G F ) 



Gp 
V2 



(96) 
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Figure 4: a) Tree- level Feynman diagram for |A5| = 1 transitions, at the lowest 
order in Gf and without strong- interaction corrections; b) the same diagram in the 
effective theory. 




Figure 5: Leading-order QCD corrections to the diagram in Fig. a) in the full 
theory; b) in the effective theory. 

As can easily be verified by expanding the W propagator in Fig. in the absence 
of QCD interactions the sum in the right-hand side of Eq. (Kfl) would be restricted 
to only two operators, 

Q± = 2 [ s L ^u L u L ^d L ± sl^lUlYul ] , (97) 

whose (Wilson) coefficients are given by 

C + = C_ = 1 . (98) 

The situation becomes more complicated when the QCD corrections in Fig. |5] are 
taken into account. In order to evaluate these effects within perturbative QCD 
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it is assumed that external quarks carry large momenta. Moreover, since some 
diagrams develop artificial divergences (absent in the real case) within the effective 
theory, it is necessary to regularize the effective operators by means of the Wilson 
coefficients, introducing an arbitrary renormalization scale [L. The requirement that 
the product Cj(/i) x Qi(fi), and thus all physical observables, be independent of /i, 
fixes the renormalization-group evolution of the coefficients Ci as a function of \x 
unambiguously. The initial value of the C, is fixed by matching the effective theory 
to the full one at high scales [i.e. for external momenta of 0(My/)\, where the 
perturbative calculation is rapidly convergent. Then, using renormalization-group 
equations the Ci are evolved down to a scale of (9(1 GeV), as close as possible to 
the physical scale of the process (the kaon mass), but still high enough to trust 
perturbative QCD. 

As a result of this procedure, all leading logarithms, i.e. all the terms of 
0[a™ \og(M w / fi) n ], are encoded in the Wilson coefficients. At this level of accuracy 
the set of effective operators is still the one in Eq. fl97|); however, their coefficients 
have changed to 0) 



c+00 




(*s(M w ) \ a "b = ( a s {M w ) \ **> 



where /3q = (33 — 2Np/12) and Np denote the number of dynamical quarks. The 
weight of the operator Q_ is therefore enhanced, whereas that of Q + is decreased. 

This program can be iterated at the next-to-leading order, resumming also 
subleading logarithms. At this level the number of terms increases, with the inclusion 
of the so-called penguin operators. Finally, when also the small electroweak 0(e 2 Gp) 
operators relevant to CP- violation studies are taken into account, the full basis of 
^e//' _1 below the charm threshold contains, in addition to Q±, the following set of 
operators: 2D 



Q3(5) 


= 4 
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q=u,d 
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Qa(6) 


= 4 




E 

q=u,d 
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Q7(9) 


= 6 




E 

q=u,d 
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^8(10) 


= 6 




E 

q=u,d 


e l ^(L)7^S(L) 

s 



(100) 
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where a and (3 are colour indices and e q denotes the electric charge of the quark q. 

Thanks to the resummation of the subleading logarithms, the coefficients of 
the Qi are known with high accuracy. ^ I^H) According to Eq. (|96|), the remaining 



problem to be addressed in order to calculate the transition amplitudes of physical 
processes is the evaluation of the hadronic matrix elements of the Qi. This issue is 
the goal of the second EQFT construction. 

4.2 Lowest-order non-leptonic |A5| = 1 chiral Lagrangian 

Following the basic principles of CHPT, in order to compute non-leptonic weak 
transitions of pseudoscalar mesons we need to construct the chiral realization of 
^e//' -1 ' or we nee d to consider the most general Lagrangian, written in terms of 
pseudo-Goldstone fields, transforming as li)^ -1 under SU(3)l x SU(3)r. 

The operators of Eq. ( |100| ) transform linearly under SU(3)i x SU(3)r in 
the following way: 

Q-, Q3, Q4, Q5, Q§ (8l, 1r) , 

Q+, Q9, Qio (8 L , lfl) + (27 L , , (101) 

Q71Q?, (8 L ,8 R ) . 

Analogously to the case of light-quark masses, chiral operators transforming like 
the Qi can be built by introducing appropriate external sources. As an example, in 
order to build the (8l, 1r) operators we introduce the source 

X^gM 1 ( 102 ) 

and we consider all possible operators invariant under G and linear in A. Then, 
fixing the source to the constant value 

A -> A = ^(A 6 -iA 7 ) = I , (103) 

we automatically select the | AS\ = 1 component of all possible (8x, 1r) terms. For 
(27 l, 1r) and (8^,, 8r) operators the procedure is very similar, the only change is the 
source structure. 

Interestingly, within each group of transformations there is only one chiral 
realization at the lowest order. These are given by: IE H) 

= (XL^) (8 L , 1 R ) 0(p 2 ) , 

= (L,) 23 (L») U + l(L„) 21 (IS) 13 (27 L , 1 K ) 0(p 2 ) , (104) 

wt ] = F 2 (XWQU) (8 L ,8 R ) 0(pP), 
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where L M = u^u^u. While the singlets under SU (3)_r are at least of 0(p 2 ), the lowest- 
order realization of (8^, 8#) operators starts at O(p ). This difference does not create 
a power-counting mismatch since at the quark level the (8^, 8r) terms are suppressed 
by a factor e 2 with respect to the SU (3)r singlets. Thus we can consistently consider 
the chiral realization of the \AS\ = 1 Lagrangian at O(GFp 2 e ) + O{Gpp e 2 ): 



G 8 W^ + G 27 W$ + GsW^ +h.c. (105) 



The constants Gi appearing in Cy) are not fixed by symmetry arguments, 
apart from the constraint ^s(Gj) = 0, which holds in the limit where CP is an exact 
symmetry. By construction we can only write 

r{I - F) = -^v: s V ud J2C^)(F\Q^)\I) = F^G t .(F\wl 2) \I) + 0(p*) 

(106) 

and unfortunately in this equation we have unknown terms on both sides of the last 
identity: the (F\Qi(p)\I) on the left and the Gi on the right. As we will show below, 
since the Gi are few, we can determine them all by comparison with experimental 
data in K — > 2ir decays, i.e. by measuring T(K — > 2tt). In this case Eq. ( |106| ) lets 
us fix the matrix element of some combination of Qi{fi) between a kaon and two 
pions and then, more interestingly, lets us predict the matrix element of the same 
combination of Qi(fi) in other channels, e.g. between a kaon and three pions. In the 
future one could hope to replace, fully or in part, the use of experimental inputs with 
a theoretical determination by means of lattice-QCD of some (F\Qi(fi)\I) . Note, 
however, that in both approaches, either using lattice QCD or experimental data, 
the role of CHPT is to relate the matrix elements of the Qi(y) in different channels 
and not to predict from scratch the (F\Qi(fi)\I) in a given channel. 

Rough theoretical estimates of the Gi can be obtained by using simpli- 
fying assumptions that let us compute the (F\Qi(ji)\I) explicitly, although with 
large (and typically uncontrolled) uncertainties. These are very useful to under- 
stand the order of magnitude of the Gi. The most natural simplifying assumption 
to estimate the (F\Qi(ji)\I) is the factorization hypothesis, that can be formally 
justified within QCD in the limit of an infinite number of colours (Nc — > oo). 
According to this hypothesis, the hadronization of operators such as Q±, with a 
colour- singlet (current) x (current) structure, is given by the product of the corre- 
sponding hadronized currents. Since the hadronization of a colour-singlet quark 
current is completely determined by the coupling of the external sources in the 
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strong Lagrangian [see Eqs. (|T^)-(^)]. 

fclyfL - ^ ( d ^\j) oi , (107) 
this hypothesis leads to a full determination of the hadronization of Q±: 



I w m + — w {2) 



(108) 



In this case there is no trace of the renormalization scale /x, and indeed it can be 
shown that the anomalous dimensions of Q± vanish in the limit Nc — > oo. Using 
Eq. ( |106| ) one then finds 

G 27 | fac , = ^ -^CK,, (109) 



where in the terms after the arrows we have employed the values of C± in Eq. ([53), 
which are consistent with the absence of anomalous dimensions for Q±. As can be 

(2} 

noted, the normalization of C)J is such that the Gi are expected to be 0{GFV us Cj), 
where Cj are Wilson coefficients of partonic operators transforming as the corre- 
sponding chiral realizations. 

As anticipated, the Gi can be experimentally determined by means of K — > 
2tt data. As long as we are interested in the real parts or the absolute values of the 
decay amplitudes we can neglect G$. Indeed the Wilson coefficients of (8l,8r) 
operators are absolutely negligible with respect to the others - for the imaginary 
parts this is not the case. Introducing the isospin amplitudes A and A 2 

A(K° -> 7T+7T-) = A e lSo + -^A 2 e iS2 , 

A(K° -> ttV) = A e iSo - V2A 2 e l52 , (110) 
A(K + -> tt+tt ) = ~A 2 e i62 , 

it is easy to verify that 

A = V2F (Ml - Ml) (g 8 + ^G 27 ) , (111) 

A 2 = 1 1fG 21 (m 2 k -MI) . (112) 
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The comparison with the experimental data then leads to: 

|G 8 | = 9.1 x KT 6 GeV" 2 , (113) 
G 27 /G s = 5.7 x 1(T 2 . (114) 

Interestingly, the experimental values of G& and G27 are substantially different from 



their naive estimates in Eq. (|109|) : the absolute value of Gg is about 8 times larger, 
whereas IG27I is reduced to 50%. This phenomenon, known as the "A/ = 1/2 rule", 
is in part explained by the running of the Wilson coefficients. Indeed if in Eq. ( |109| ) 
one uses the leading-log values of C±(fi), with a renormalization scale fi ~ 1 GeV, 
then the factorized estimate of \G$\ increases by a factor of about 2 and the one of 
IG27I reduces to 70%. Although this is encouraging, it is clear that there still is a 
large non-perturbative effect hidden in the matrix elements of four-quark operators, 
especially the (8l, Ir) ones. At the moment the best we can do is to measure this 
effect using K — > 2rr data. 

Once the Gi have been fixed from K — > 2tt, the theory is absolutely pre- 
dictive in all the other non-leptonic channels. In Table |2| we show the comparison 
between theory and experiments in K — > 3n amplitudes. The latter are classified 
according to the variation of isospin (AI = 1/2 or AJ = 3/2) and to the depen- 
dence from Dalitz-plot variables, as summarized in Table [| In agreement with naive 
power counting, the discrepancy between lowest-order chiral predictions and data, 
within the dominant amplitudes, turns out to be around 30%. On the other hand, 
0{p 2 ) operators have not enough derivatives to produce a non-vanishing result for 
the suppressed quadratic slopes. The comparison between data and lowest-order 
CHPT predictions is much more satisfactory in the case of K — > decays, where 
the 0(p 2 ) amplitudes coincide with the bremsstrahlung from K — > 2tt, which is 
known to be largely dominant. HH) 

4.3 0{p A ) counterterms 

In order to obtain a more refined description of experimental data, able to include 
absorptive effects and subleading amplitudes (such as K — > 3tt quadratic slopes and 
K — > 7T7T7 direct-emission terms) , it is necessary to go beyond the lowest order in the 
chiral expansion. Since we are interested only in contributions of order Gf, we can 
proceed analogously to the case of strong interactions with the simple substitution 

r (2) _ r (2) r (2) 
*-S S '-W 1 

4" - 4" + 4?. (n5) 



3S 



Amplitude 


o(p 2 ) 


G(p 4 ) 


Experimental fit 


«i 


74.0 


input 


91.71 ±0.32 


/3i 


-16.5 


input 


-25.68 ±0.27 


Ci 





-0.47 ±0.18 


-0.47 ±0.15 






-1.58 ±0.19 


-1.51 ±0.30 




-4.1 


input 


-7.36 ±0.47 




-1.0 


input 


-2.43 ±0.41 


73 


1.8 


input 


2.26 ±0.23 






-0.011 ±0.006 


-0.21 ±0.08 






0.092 ±0.030 


-0.12 ±0.17 






-0.033 ±0.077 


-0.21 ±0.51 



Table 2: Experimental results vs. CHPT predictions in K — > 3ir decays. EUP 





AI — 1/2 


A/ = 3/2 


Phases 


Constant terms 


1 («i) 


i M 


1 


Linear slopes 


1(A) 


2 (Aj,7 3 ) 


3 


Quadratic slopes 


2 (£,&) 


3 (Cs,6,&) 




0(p 4 ) free parameters 


2 


3 





Table 3: Number of independent isospin amplitudes vs. 0(p A ) free parameters in 
K — > 37r decays. 



Here denotes the most general C(p 4 ) Lagrangian transforming like Cy) under 
chiral rotations and thus able to absorb all the one-loop divergences generated by 

The C(p 4 ) operators transforming like (8l, 1_r) and (27/,, 1r) were classified 
for the first time by Kambor et al, 23) about 10 years ago, whereas the 0(e 2 p 2 ) 
terms transforming as (8/, 8#) have been analysed only very recently. HH) The overall 
picture is certainly worse than in the case of strong interactions, since the number of 
independent operators is much larger: already within the dominant (8^, 1r) sector 
there are 37 independent terms. HIP Nonetheless the theory still has a considerable 
predictive power, since out of these terms only few combinations appear in observable 



processes. B 0) For instance, within K — >• 3ir decays the number of independent 



counterterm combinations coincides with the number of leading (constant and linear) 
amplitudes. Thus a fit of the C(p 4 ) parameters from the leading amplitudes leads 
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Figure 6: 0(p A ) diagrams for the transition Ks — > 77. The black box denotes the 
weak vertex. 

to unambiguous predictions for the quadratic slopes, as shown in Tables ^| and |3|. 

There are even cases of physical processes that receive 0(p 4 ) contributions 
only from loop diagrams and not from counterterms, like the Ks — > 77 decay The 
leading contributions to A(K S — > 77) are the loop diagrams in Fig. ||. Since there 
are no local contributions to this transition at 0(p 4 ), the loop amplitude turns out 



to be finite 23) and leads to the following parameter-free result: 

B{K S -> 77) (4) = 2.1 x KT 6 . (116) 

The above prediction turns out to be in excellent agreement with the experimental 
data, B(Ks — > 77) exp - = (2.5 ± 0.5) x 10~ 6 , providing a very significant test of the 
quantum nature of this effective field theory. 

4.4 Beyond 0(p 4 ) 

Another process that receives 0(p 4 ) contributions only from loop diagrams is Kl — > 
7r°77. The diagrams describing this transition are very similar to those in Fig. |6|; 
the only difference is an external pion leg attached to the weak vertex. In this case, 
however, the 0(p 4 ) parameter-free prediction of the branching ratio 

B{K L -> A7) (4) = 0.6 x 10~ 6 (117) 
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turns out to be in bad agreement with the experimental finding: B(Kl — > 77) cxp - = 
(1.7 ± 0.1) x 1CT 6 . What is the reason for the big difference between K$ — > 77 and 
Kl — > vr°77? The answer to this question can be traced back to what we learned in 
lecture 3. 

The first difference is related to the absorptive parts of the amplitudes. 
Since the coupling G$ has been fitted from K$ — > 2tt, the imaginary part in 
Kg — > 77, computed at the one-loop level, perfectly reproduces the absorptive 
contribution due to the Ks — > 2tt intermediate state. On the contrary because 
of the underestimate of K — > 3n amplitudes at 0{p 2 ) (see Table §), the one-loop 
imaginary part of Kl — * 7r°77 underestimates the real absorptive contribution to 
this channel by 20%-30%. The second important difference between K$ —>■ 77 and 
Kl — > 7r°77 is induced by resonance contributions. Vector and axial-vector mesons, 
which are known to produce sizeable effects in the strong sector (see Table [D, do not 
appear in Kg — > 77 but can affect the Kl — > 7r°77 amplitude. In the latter case one 
can therefore expect a sizeable local 0(p 6 ) counterterm, encoding the contribution 
induced by vector-meson exchange. 

These two effects shows that, contrary to Ks — > 77, a good description of 
Kl — > vr°77 requires the inclusion of 0(p 6 ) terms. 0) A similar situation emerges in 
many other channels and can be interpreted as a general rule: unitarity corrections 
induced by pion loops and vector meson resonances provide a useful guide toward 
the estimate of the most significant higher-order corrections. 

5 Conclusions 

Solving strong interactions at low energy is a very difficult task, which so far nobody 
has been able to complete with analytical methods. On the other hand, at low 
energy, the approximate chiral symmetry does impose severe constraints on Green 
functions in the form of Ward identities. In studying the phenomenology of strong 
interactions at low energies it is extremely useful to take into account these symmetry 
constraints. The effective Lagrangian method is a tool to derive these constraints 
in an automatic manner - at the same time respecting also the general properties 
of unitarity and analyticity. 

In this series of lectures we have introduced the basic concepts and technical 
tools for using this method. The last lecture was devoted to the study of a few 
examples of physically interesting decay channels of the kaons. In recent years a rich 
experimental activity at various facilities around the world has started, with the aim 
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to study, with very high accuracy, all possible decay modes of the kaons, in some 
cases measuring extremely small branching ratios. The use of chiral perturbation 
theory for analysing this rich phenomenology is essential: it is not by chance that 
this lecture series was given at the Frascati Laboratories, which host one of the main 
world facilities for this kind of physics. 

Kaon physics is certainly not the only subfield of strong interactions where 
CHPT is useful and is applied successfully. A subject that we have not touched 
at all in these lectures, but where the application of CHPT is nowadays routine 
is the physics of baryons, going from the "easy" one-baryon sector (irN scattering, 
ttN — > inrN, photo- and electro-production reactions, etc.) to the more complicated 
sectors involving two or more baryons. 

In perspective, we believe that CHPT will become more and more used 
in connection with lattice calculations: as we have seen in the few examples of 
calculations discussed here, the dependence on the quark masses is always explicit. 
Whenever a quantity is calculated on the lattice one can use CHPT to determine the 
explicit dependence on the quark masses and make the extrapolation to the physical 
quark masses in a controlled manner. The "simple procedure" we described here 
hides in fact some serious technical difficulties - conceptually, however, it is obvious 
that merging the two methods offers clear advantages, and we will see more and 
more of them in the future. The technical difficulties are only welcome, because it 
means there will be work for theorists! 

Let us finally conclude by mentioning that this method is also applicable 
to very different physical systems. In fact it can be applied to any physical system 
in which spontaneous symmetry breaking occurs, and which has an energy gap 
such that, at low energy, only the Goldstone modes dominate the physics. We find 
very interesting, for example, the applications to magnetic systems - in some cases, 
precisely the same effective Lagrangian that we have constructed for pions can be 
used to describe the behaviour of one of these. H§D 
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